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Abstract 

In this paper, incremental exponential asymptotic stability of a class of switched Caratheodory nonlinear systems is studied based on 
the novel concept of measure of switched matrices via multiple norms and the transaction coefficients between these norms. This model 
is rather general and includes the case of staircase switching signals as a special case. Sufficient conditions are derived for incremental 
stability allowing for the system to be incrementally exponentially asymptotically stable even if some of its modes are unstable in some 
time periods. Numerical examples on switched linear systems switching periodically and on the synchronization of switched networks of 
nonlinear systems are used to illustrate the theoretical results. 
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1 Introduction 


Studying incremental stability of nonlinear systems is par¬ 
ticularly important in many application areas, including 
observer design and, more recently, consensus and synchro¬ 
nisation problems in network co ntrol where conv e rgence 


analysis is a fundamental sten (IWans 8c Slotinel 

2005 

Russo & di Bernardol l2009allbl iRusso et al.. 20101 

2011 

201,31). 



Since the early work bv iLewisI (ll949h : lDemidovi^(Il967h . 
contraction theory has been highlighted as a promising ap¬ 
proach to study incremental exp onential asymptotic s t ability 
((5EAS) of nonlinear sys t ems (iLohmiller & Slotinel 1993 
iForni & Sepulchr^ 120141; lAngelil l2002h : also see JouffrovI 
( 2005ll for an historical ove rview. In particular, as shown by 
ILohmiller & Slotin3 (1 19981) . sufficient conditions for (5EAS 
of a given nonlinear system over an invariant set of interest 


can be obtained by studying the matrix measure of its Jaco¬ 
bian induc ed by some vector norm. It is possible to prove , 
as done bv lLohmiller & Slotinel (I1998I) : IRusso et aD(l2010h . 
that if such measure is negative definite in that set for all 
time then any two trajectories will exponentially converge 
towards each other; the rate of convergence being estimated 
by the negative upper bound on the Jacobian measure. 

Numerous applications of contraction analysis have been 
presented in the literature from observe r design to the synthe¬ 
sis of network control systems. Se e e.g.ILohm iller & Slotinel 
(119981) : iForni & Senulchr^ (I2014|) :Ir usso et a IJ (120111) . Re- 
markably, the problem of studying incremental stability of 
switched and hybrid systems has attracted relatively lit¬ 
tle attention despite the large number of potential appli¬ 
cations, e.g. power electronic networks, variable structure 
systems, walking and hopping robots, to name just a few 
ddi Bernardo et al.L 120081 : iLiberzoni l2003t ICortesL l2008l) . 
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It has been suggested by iRusso & di Bernardol (1201 ll) : 
Idi Bernardo et al. (12014l) that extending contraction anal¬ 
ysis to this class of systems can be a viable and ef¬ 
fective approach to obtain conditions for their incre¬ 
mental asymptotic stability. Related approaches include 
the work on convergence of piecewis e affine c ontinu ¬ 
ous systems presented bv iPavlov et alJ (l2005bllaL l2007l) : 
iPavlov & van de Wouwl ( 20081) and the recent conference 
papers ddi Bernardo & Liuzzal 120131 : Idi Bernardo & Eio^ 
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12014 . 

One limitation of the existing e xtensions of contraction the- 
ory to switched systems, (e.g. Idi Bernardo et al.L I20l4 . is 
that they rely on the use of a unique matrix measure to as¬ 
sess the Jacobian of each system modes. This is a particu¬ 
larly restricting assumption as it would be desirable to use 
measures induced by different norms to evaluate the Jaco¬ 
bian of each of the system modes. This would correspond to 
studying incremental stability of the switched system with 
multiple incremental Lyapunov functions rather than using 
a common one (which is much harder to find). 

The aim of this paper is to address this problem and present 
conditions for contraction and incremental stability of a large 
class of switched Caratheodory systems. The key idea is 
to define a novel concept of matrix measure via multiple 
norms and exploit the transaction coefficients between these 
norms. In so doing, sufficient conditions are derived for in¬ 
cremental stability that allow for a system to be bEAS even 
if some of its modes are unstable (or not contracting) over 
some time intervals. The theoretical results are illustrated 
via their applications to some representative examples, in¬ 
cluding synchronisation in blinking networks. 


2 Preliminaries 

We focus on switched dynamical systems of the form 

x = f{x,r{t)) (1) 

where x G M" and the switching signal r{t) is assumed to be 
a real-valued piecewise continuous (PWC for short) function 
with respect to time: there exist countable discontinuous 
points to < ti <■■■< ti <■■ ■ such that r(L±) exist and 
r{ti) = r{ti+) for all ti. A typical example is the staircase 
function r{t) = for ti < t < L+i, * = 0, 1 , • • •, for the 
increasing time sequence {tj}j>o, which has been widel y 
used as switching signal in control systems (iLiberzoni l2003h . 

Here we make the following hypothesis: 

TLi'. For some C C R", the vector field f{x,r{t)) : C x 
R>o —>■ R" is (i) continuous with respect to (a;,r); (ii) 
continuously differentiable with respect to x, and (iii) 
there exists a Lebesgue measurable function m{t) such 
that \f{x, r(f))| < m{t) for all a; € C and t € R>o- 

It can be seen that under hypothesis TLi, and with r{t) 
defined as above, the vector field f(x,r(t)) defines a 
Caratheodory switched system (iFilippovl Il988h . It can be 
proved that, given an initial condition in C, a solution of a 
Carathedory system exists and is unique (iHalel 11954 . We 
define (pit] to, Xq, rt) as the solution of (1) with x{to) = xq 
and the switching signal r(f), where rt denotes the trajec¬ 
tory of r{t) up to t, i.e., r* = {r(s)}t(,<s<t. 


In this paper, | • 1^ stands for a specific vector norm in 
Euclidean space and the matrix norm induced by it, which 
can be defined in different ways that are all equivalent. The 
transaction coefficients f rom the norm | ■ I „ to | • |b is defined 
as f3ab = |x|t, jBourbakii 119781) . 

A continuous function a : [0, a) [0, oo) is said to belong 
to class K. if (I) it is strictly increasing; (II) q:(0) = 0. And, 
a continuous function ffip^t) : [0, a) x [0, oo) —>■ [0, oo) is 
said to belong to class ICC if (1) for each fixed t, the function 
/3(p, t) belongs to class /C; (2) for each fixed p, the function 
Pip, t) is decreasing with respect to t and limt^oo Pip, t) = 
0. In addition, if a function Pip,t) of class ICC converges 
to 0 exponentially as f — oo, Pip,t) is said to be of class 
8ICC. Here, we give the foll owing definition of incremental 
stability from lAngelil (l2002h with modifications. 

Definition 1 System (1) is said to be incrementally asymp¬ 
totically stable (6 AS for short) with r(f) in the region C C 
R" if there exists a function Pis,f) of class ICC such that 
for any initial data Xo,yo G C and starting time to, the fol¬ 
lowing property holds 

\(j)it + to;to,xo,rt) - (pit -G to;to,yo,rt)\ < Pi\xo - yo\,t) 

for some norm \ ■ \. If Pis, t) is picked independently of the 
initial time to, then system (1) is said to be incrementally 
uniformly asymptotically stable iSU AS for short). If Pis, t) 
is of class 8ICC, then system (1) is said to be incrementally 
exponentially asymptotically stable (SEAS for short) and 
incrementally uniformly exponentially asymptotically stable 
(SUE AS for short) if Pis,t) is chosen independently of to- 

Definition 2 A set C C R" is said to be a n-reachable 
set if there exists a continuously differentiable curve 7(5) : 
[0, 1 ] —>■ C that links xq and yo, i.e., 7 ( 0 ) = Xq and 7 ( 1 ) = 
yo, and satisfies |7'(s)lx(to) < - yolxM ’® ^ 

[0, 1 ] and some constant k > 0, independently of Xo and yo. 

3 Switched matrix measures and general contraction 
analysis 

The matrix measure induced by the vector norm | • where 
X is the index for the norm being used, is defined as 

= lim T [\iln + hA)\^ - 1] 

n^0+ ri 

for a square matrix A G R"’” and was used for the con- 
traction analysis of smooth nonlinear systems, see e.g. 
lLohmiller& Winel(l1998h . 

Given a PWC function x(f), the left (right) limit of | • at 
time t is defined as lim/i^o- \x\x{t-vh) (lim^^o-r \x\xit-vh)) 
if it exists for all x G i?", and is denoted by | • \x{t±) respec¬ 
tively. We say that the switched norm | • |^(t) is continuous 
at time t if |a;|j^(t_) = |a;|^(t+) = |a;|^(t) for all x G R”, i.e. 
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if I • |^(t) is left and right continuous. We say that | • is 
uniformly equivalent if there exists a constant D > 0 such 
that < D\x\y^(^s) for x G R" and f, s G R. 

We can now extend the definition of matrix measure to the 
case of multiple norms, taking the time-varying nature of 
x{t) into consideration, as follows. 

Definition 3 The switched matrix measure with respect to 
multiple norms \ ■ is defined as 

^x(i)(^) = 1 ™ T ^ hA)x\^(^t+h) - ^] ( 2 ) 

h^o+h = 1 


if the limit exists, where lim stands for the limit superior. 

It can be seen that if x{t) is constant over an interval, say 
[f, t + 6), then (A) = (^) over that interval. 

The existence of the switched matrix measure is related to 
the partial differential of the switched norm | • defined 
as follows 


9t{\ ■ lx(t)) 


lim sup 


\^\xit+h) - 1 

h 


We say that the multiple norm | • |^(t) is right regular at time 
f if i9t(| • lx(i)) exists. Thus, we have 

Proposition 1 If\-\j.^f-^is right regular at t, then (i). | • | 
is right continuous at t; (ii). exists at t. 


Lemma 1 Suppose that \ ■ |^(t) is right-regular in 
Consider the following time-varying Caratheodory dynam¬ 
ical system x{t) = A{t)x(t) with some PWC matrix-valued 
function A{t) G R"’” and x{t) G R". If < aft) 

for some measurable function aft) for f G [ti, t 2 ], then 


\x{t2)\x(t2) < exp 





(3) 


Proof For any t G (^ 1 ,^ 2 ) except discontinuous points of 
Aft), consider the following quotient 


I{h) = ^[\x{t + h)\^(^t+h) - \x{t)\xit)] 

= Ms)xfs)ds\x(t+h) - |a;(f)lx(i)] 

=+ hA(t))xft) [^(s)a;(s) 

-A{t)xft)]ds\^^t+h) - |a:(f)lx(i)} 

< + hA{t))x{t)\^(^t+h) - |a;(f)lx(t) + o{h)], 


which implies that lim;i^of^(^) < ’^x(t)iMt))\x{t)\x{t)- 

Thatis, D+\x{t)\^(^t) < ’^xit)iMt))\xit)\x{t) < a(f)k(f)lx(t) 

holds for almost every t G [ti, ^ 2 ], where D+ stands for the 
Dini derivative. Thus, Ineq. (3) holds, noting that xft) is 
continuous with respect to t. ■ 


We make the following hypothesis on the multiple norm. 


Proof The first statement is straightforward from the defi¬ 
nitions of 9t(| • lx(t)) ^nd right continuity of | • lx(t)- 


T-L 2 ' I • lx(t) is right-regular everywhere but at the time instants 
{tj}, it is right-continuous and uniformly equivalent, and 
its left-limit exists at each h. 


The quotient term of the definition of iy.^^^t^{A) gives 


< 


- sup [|(/„ -b hA)x\^(^t+h) - \{In + hA)x\^(^t) 
^ l^lx(t) = l 

+ |(Fn + h.A)x\^(^t) ~ 1 ] 

Stip [|(-^n “b h'A)x\^^t^ii'j |(fn “b ^^)^lxC*)] 
l^lx(t) = l 

+ T sup [|{/„ -b hA)x\^(^t) - 1] 


- f ■ lx(t))K-^" + ^^)lx(t) + + o{h)] 

where o{h) is an infinitesimal term with limh^o o{h)/h = 
0. So, Vx{t)i^) ^ Fx(t)(^) +^t(l • lx(t)) (statement (ii)). ■ 


Thus, we denote by N{t, s) = ff{i : s < ti < t}, the num¬ 
ber of time instants ii in the time interval (s,f); obviously, 
Nft) = N{t, to). Here, ff stands for the cardinality of a set. 
In addition, we make the following assumption. 

"Ha: There exists a K-reachable set (IRusso et al.L 1201 Oh C C 
R" which is a forward-invariant for (1). 

Then, we have a general result on SEAS of switched dynam¬ 
ical system (1) by contraction analysis in multiple norms. 

Theorem 1 Suppose that hypotheses "Hi, 2,3 hold and that 
rft) is PWC. If there exist a measurable function a(t), non¬ 
negative constants fSj, j = 1,2, •••, c > 0, and Tq > 0, 
such that, for all x € C, the following conditions hold 

Vx{t) <oi.{t),\Itf^ij, x&C 

I ■ Ixfe) Ixfe-)’ ^ d 


Using the definitio n of switched m atrix measure, we extend 
Coppel inequality (ICoppelill96-5h as follows. 
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(4) 

(5) 












and, for all T > Tq. 


we have 



then system (1) is SEAS in C with respect to r(t); if (6) 
holds for some c > 0 independent of then (1) is SEU AS 
with respect to rif) in C. In addition, the exponential con¬ 
vergence rate can be estimated as 0 (exp(—c(t — to)))- 


|w(t,s)|^(t) < exp 
< exp 


a{T)dT 




/■* 


/ a(r)dr 

.N{t) 

■ exp 


^N{t) 


'to 


a{T)dT |w(to,s)|x(t„). 


Proof For any Xo,yo € C, C being K-reachable for some 
K > 0 implies, from Definition 2, that there exists a con¬ 
tinuously differentiable curve 7 ( 3 ) : [ 0 , 1 ] —?► C such that 
7 ( 0 ) = xo and 7 ( 1 ) = yo, and | 7 '(s)lx(io) ^ - 

yolx(to)’ for ^11 ■s ^ !]■ 


Let = (j){t-,to,y{s),rt), s G [0,1], be the solution 

of (1) with initial value 'ijj{to,s) = 7 ( 5 ). Since f{x,r{t)) 
is continuous with respect to {x, t) except for the switching 
time points {tj} and continuously differentiable with respect 
to X, then to,xo,rt) is continuously differentiable with 
respect to the initial value xq. Let xq = 7 ( 5 ). Then, w = 
^ is well defined and continu ous. By the same algebraic 
manipulations first presented by iRusso et"^ (1201 Oh . except 
for {tj}, w{t, s) is the Caratheodory solution of: 

|w;(fo,s) = 7'(s) 


Consider the sequence \w{tj, s)|^(t .) of Eq. (7). For all s G 
[0,1], by the extended Coppel inequality (3) and condition 
(4), we have 


< exp 


i{t)dt 


' tj —I 




Using inequality (5) between | • and | • gives 


< Pj exp 


a{t)dt 


k(fVi>s)lx(t,_i) 


Then, iterating down to z < j, we have 


\w 


^feexp 


J a{t)dt 


k(f*:S)lx(t7- 


For any t, noting that from definition fjv(t)+i > t > fjvftp 


‘'Ar(t)’ 


Hence, we obtain 


\w{f s)lx(t) < exp(-c(f - fo))|7'(s)lx(to)> ( 8 ) 

for all t > to Tq. Since c is independent of s, SEAS 
ca n be derived by foll owing similar arguments as presented 
bv iRusso et al.l (l2010h . In detail, since | • |^(t) is uniformly 
equivalent, there exists D > 0 such that |a;|^(t) < 
for any x G M" and t, t' > tg. Then, by ( 8 ), we have 


lf(t;x(0),rt) - f(t;y(0),rt)l^(to) 
< ^l(/>(t; x(0), rt) - fit; y(0), rt)|x(t) 
diljit,s) 


(9) 


= D 


ds 


-ds 


< D 


x(t) 


wit,s)\^i^t)ds 


<D f exp(-c(f-fo))l7'(s)lx(io)r^s 
Jo 

= Dnexpi-cit - to))\xo - t/oIxOo) 


( 10 ) 


which converges to zero exponentially. This proves that sys¬ 
tem (1) is SEAS. 


If ( 6 ) holds for some c > 0 independently of to, then ( 8 ) 
holds independently of to if t — to > To, which implies 
limt_>oo l'tf(0lx(<o) = 0 is uniform with respect to to- By 
the arguments above, system (1) is SUE AS. In addition, 
inequality ( 10 ) implies that the exponential convergence rate 
can be estimated as 0 (exp(—c(f — fo)))- ■ 


Note that c is an estimate of the exponential convergence 
rate and also an index of the average contraction rate of ( 1 ). 

As a simple example to illustrate this result, define a 
switched matrix measure based on the following multi¬ 
ple norms: |x|p(t) = y/x^P{t)x. Here, P{t) G K”’" 
is a differentiable matrix-value function where each 
P{t) is a symmetric positive definite matrix such that 
its largest eigenvalue Amax(-P(f)) the smallest eigen¬ 
value XjniniPit)) is upper bounded and lower-bounded 
positive respectively. Then, for any A G R"’", by 
simple algebraic manipulations, it can be shown that 
the induced switched matrix measure is z/p(t)(A) = 

^An.ax[P-l/2(f)(-P(f)4l + A^P{t) + P(<))P-l/2(t)]. 
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Now, consider the linear time-varying (LTV) system x = 
A{t)x + I{t) with some PWC matrix function A{t). Then, 
(6) in Theorem 1 can be fulfilled by assuming that it holds 
that P{t)A(t) + A(t)^P{t) -f P{t) < 0 for almost every t. 
A similar condition can also be obtained by using the theory 
of convergent systems (iPavlov et all 1200% . 


Remark 1 Analogous arg ume nts to those prese nted by 
\Lohmiller & SlotinA il99a) and \Russo et al\ ( 120701) can be 
followed to prove that the contracting system converges 
towards a periodic solution if r(t) is periodic. Indeed, in 
this case the time instants where rit) switches, say {tj}, 
are also periodic with period T and so is the function 
N{t, s) = N{t + T, s + T) for all t > s. Then, inequality 
(6) can be fulfilled by assuming that 


1 

T 



^log(/3j ) 

7=1 


< —c. 


( 11 ) 


where J is the number of switches in one period. (Note that 
this inequality holds independently of the initial time to, due 
to the periodicity.) From Theorem 1, one can conclude that 
system (1) is SEAS. The existence of a periodic solution of 
(1) follows by using the contraction mapping theorem. 


4 Contraction analysis of staircase switching and trans¬ 
action coefficients between norms 


As a special case and an important application, in this sec¬ 
tion, we assume that the switching signal r{ t) of the switched 
system (1) is a staircase function: 

Pi: r{t) = ii, ti < t < L+i, i = 0,1 ,---, with an 
increasing time point sequence {ti} with ti = 0 and 
limi_^oo ti = oo, where € S where S is a countable set. 


Let Ai = ti+i — ti > 0 he the interval between two points, 
and N{t, s) = : s < ti < t}, be the number of time 

instants when r(t) switches that fall in the time interval [s, t], 
in particular, N(t) = N{t,to). 


By using multiple norms | • |^(t) with x{t) = r{t), we obtain 
the following result as a consequence of Theorem 1. 


Corollary 1 Suppose that hypotheses 771,2,3,4 hold. If there 
exists constants ai, nonnegative constants jdi, i = 1, 2, • • •, 
c > 0 and Tq > 0 such that, for all i and x € C, the 
following conditions hold 

Tii 0)) < aj, I • <Pj\- , V j (12) 


and, for all T > Tq, 



then (1) is SEAS with respect to r(t) in C; if (13) holds for 
some c > 0 independently of Iq, then (1) is SUE AS with 
respect to r(t) in C. Moreover, the exponential convergence 
rate is estimated as 0(exp(—c(f — fo)))- 


Proof Let a{t) = ai if t G We have 




Note that 

rT 

'to 


1 1 

tL ? 


— f^r{t) 



< a{t), M tfif. 


O^N{to+T) + l(to + T — fiv(to-l-T)) + 

N{to+T) , 


+ cviAi 


i=0 

Under condition (13), one can conclude that there exists 
Ti > 0 such that 


1 

T 



N(to+T) 

a(t)dt + ^2 A 

i=0 


< —c 


(14) 


for all T > Ti. By employing Theorem 1, SEAS of (1) 
can be proved. Furthermore, if (13) holds independently of 
the initial time Iq, and therefore (14) holds independently 
of Iq, then (1) is SUEAS with respect to r(f) in C. The 
exponential convergence rate can be derived as done in the 
proof of Theorem 1. ■ 


We can exploit Corollary 1 when considering the case where 
r{t) (^i) takes values in a finite set O = {1, • • • ,K}. Specif¬ 
ically, define for any 7 > s > 0 

Tk{s, 7) = {r e [s, 7] : r(T) = fc}, 

Afki{s,t) =#{i : r{U-) = fc, and r{U+) =1, s <U < 7}, 

and constants ak and (3ki > 0, k,l = 1, ■ ■ ■ ,K such that 

Mfe <ak, I • |( < Pki\ ■ \k- (15) 

Then letting a(7) = ai when 7 G [7i,7i+i), it can be seen 
that inequality (13) in Corollary 1 can be derived if the 
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following condition holds for all T > Tq: 


for some c > 0. 


1 

T 

< 


K r K 

ctfcTfc {to, to + T) + ^ log(/?fc;)A4/ {to, to + T) 


k=l 

—c. 


1=1 


(16) 


Alternatively, condition (16) can be replaced by the assump¬ 
tion that there exists some Ti > 0 such that for each Ti- 
length interval, [nTi, {n + l)Ii], it holds that 


We now discuss the relationship between different classes 
of vector norms in greater detail. 

Quadratic norms. The quadratic norm is defined as |a;|p = 
\/Px for some positive definite matrix P. For another 
|x|q = x^Qx with symmetric positive definite matrix Q. 
Then, we have 


1 ^ 
Ti ^ 


k=l 


cikTk {nTi, {n + l)Ti) + 


K 


log(/3fei)A4i {nTi,{n + l)Ti) 


1=1 


< -c, Vn > 0. (17) 


I^Iq < v/Vax(P-l/2Qp-l/2)|:E|p (20) 

for all X G K”, where Amax(^) denotes the largest eigen¬ 
value of a symmetric square matrix A. The proof of this 
statement is straightforward and omitted here for the sake 
of brevity. 


Remark 2 From Ineq. (17), in the case that all Uk < 0, 
k = 1, • • ■ , K, which implies that all subsystems are con¬ 
tracting and hence incrementally stable, the switched system 
is incrementally stable if the duration of the mode in which 
each subsystem is active is sufficiently long. 

More specifically, if all norms | • are identical - and simply 
denoted by | • | - then Pki = 1, and inequality (13) can be 
fulfilled by simply assuming that akT{t,t -f 

T) < — c holds for all f > 0. Thus, one can see that if offc < 0 
holds for all fc = 1, • • • ,K (with respect to the same norm), 
then (1) i s incrementally stable , whic h coincides with the 
results bv iRusso & di Bernardol (1201 ih . 


A a first example, we take A(l) = 
/O -11^ 


'0 -r 

.2 -3, 


h A{2) = 


,2 -33; 


I. It is easy to see that, using the matrix mea¬ 


sures tJ-i{-) induced by the quadratic weighted vector norm 
l^^loi •= \Qix\ 2 , * = 1,2 respectively, the matrix measures 
Pi [A{i)] with 


01 


/0.707 0 . 447 ' 
V0.707 0.894 


-1 


02 


/ 0.998 0 . 322 ' 
VO.O 6 I 8 0.947 


-1 


The conditions of Theorem 1 and Corollary 1 as well as (16) 
depend on two quantities; the matrix measures ak of the 
Jacobian of the dynamical system for each constant value 
of r{t), and the transaction coefficients jdki between norms. 
As system (1) is defined in a finite dimensional Euclidean 
space, these vector norms are equivalent, i.e., there exist 
positive constants Pki such that |a;|i < flki\x\k holds for all 
X € K". The role of such transaction coefficients will be 
further investigated below. 


are negative for i = 1 and * = 2. In particular, pi [A(l)] < 
— 1 := a\ and p 2 [A{t, 2)] < —0.6807 := a 2 . 

From the results reported in Sec. 2, we get from Ineq. (20) 
that: |a:| 0 i < 1.796|a:|e2 ^^d |a:|e 2 < 1.05|a:|ej; there¬ 
fore condition (5) remains satisfied with P 12 := 1.796 and 
(321 ■= 1.05. Assuming a switching frequency of ipr = 1 Hz 
between the two modes, we have A; = 1 s for all i. There¬ 
fore (19) is satisfied and (18) is incrementally stable in 


For illustration, consider the FTV system 

x{t) = A{r{f))x{f)-\-B (18) 


where r{t) : [0, c»[h^ C = {1, 2} with A = A{1) or A{2) 
and iJ G We assume that the linear system is switched 
between two constant matrices periodically with identical 
frequency pr- Consider two vector norms | • |i. 2 , correspond¬ 
ing to matrix measures pi .2 respectively, with transaction 
coefficients / 3 i 2 and ^ 21 - Corollary 1, specifically condition 
(17), yields that the following inequality is a sufficient con¬ 
ditions for (5UEAS of system (18): 


\ [Fi{M^)) + MM‘2)) + ‘Pr ■ (log^i2 +log/32i)] 

< —c 


The evolution of the norm of the error for two trajectories 
starting at initial conditions [0.5,0.1] and [0.4, 0.1] respec¬ 
tively is shown in Fig. 1 when = [0 0] and i? = [1 1] re¬ 
spectively. In both cases we observe incremental stability as 
expected. Note that the use of two different weighted norms 
makes proving incremental stability much simpler than if 
one common metric had to be used for both modes as re¬ 
quired by previous result s (e.g. lRusso & di Bernardol 1201 It 
Idi Bernardo et all 12014h . 


/-1.3481 -2.9306\ 

Next, we take A{1) = , A{2) = 

\-2.4538 -1.2755/ 

/-11.2237 7.0628\ 

.It can be seen that A{1) is Huwitz 

\-1.7413 1.5119/ 


(19) 
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Fig. 1. Evolution of the Euclidean norm of the error \\x{t) — y{t) || 
when B = [0 0]^ (top panel) and B = [1 1]^ (bottom panel). 



Fig. 2. Dynamics of the average error Err{t) over M = 10 
independent realizations of random initial values that are 
picked from [—10,10]^ following the uniform distribu¬ 
tion. Err{t) = i - xiY + [xl - with 

Xu = {l/M) "'ith u = 1,2, where x\ 2 stands for the 

state components of the m-th realization. 

stable i.e., with all eigenvalues having negative real parts, 
but A{2) is unstable, i.e., with some eigenvalues possessing 
positive real parts. Choosing ©12 as follows: 

/0.3797 0.006l\ / 0.0644 -0.1475\ 

01 = h ©2 = 

VO.OO6I 0.4534/ V-0.1475 0.8267 / 

we find that the matrix measures, /ii(A(l)) and ^ 2 {A{ 2 )), 
induced by the 2-norms | • [e^ and | • [02 are: ^i{A{l)) = 
—2.6178, ^2{A{2)) = 0.9188. The transaction coefficients 
between the two norms | • [gj 2 can be calculated by (20) as 
I • lei < /3 i 2| • 102 and | • |e2 < /32i| • |ei with P 12 = 1.9079 
and/321 = 10.4207. The LTV system switches between these 
two modes with an identical frequency (pr (jpr = 0.25 Hz in 
this example). Then (19) holds for c = 1.1010. Therefore, 
the LTV system (18) is 5UEAS, as shown in Fig. 2, despite 
one of its modes being unstable. 


the results presented by iPorfiri et al.l (l2008h cannot be ap¬ 
plied for this specific situation. Nevertheless our extension 
of contraction analysis to switched systems gives a simple 
and viable set of conditions that can be used to prove that 
indeed the system is incrementally stable. 

More specifically, from the properties of matrix measures, 

l/2{pi{A{l)) + p{A{2)))> p{A^) 

> max{77.] (A) : A G a{{Am))} > 0 

it is not possible to find a uniform norm such that the average 
of the matrix measures of the switched matrices induced by 
this matrix norm is negative as required by condition (19). 
Therefore, in this case, multiple norms must be utilised for 
proving contraction and AEAS of the system. 

Weighted Lp-norms. The weighted Lp-type norms with 
1 < p < 00 are defined as follows. 

• Weighted Lp-norm: |a;|j,p = for some 

p > 1 and / = [Cij • • ■ ) with > 0 for all i = 
I,-- - ,n; 

• Weighted Loo-norm: |a:|^_oo = maxi^/a;/ for some ^ = 
[Cd • • • , with > 0 for alH = 1, • • • , n. 

Their transaction coefficients are summarised by the follow¬ 
ing proposition. 

Proposition 2 For p > <? > 1 with p possibly equal to 
00 and two component-wise vectors ^ = [■Ci, • • • cind 

rj = [pi, • • • , with > 0/or all i = 1, ■■■ ,n, the 

£l/P 

following hold: (1) |a;|p_5 < max^ (2) \x\q^ri < 

max^ - n^^P)\x\p^^; (3) |a:|2.5 = I^Ih with E = 

diag[Ci,--- ,Cn]- 

This result can be directly derived from iBourbakil (Il978h . 
Combining 20 and Proposition 2, we can derive all trans¬ 
action coefficients between all quadratic norms, | • |q with 
positive definite matrix Q, and all | • p for +00 > p > 1. 
The same transaction coefficients hold for the equivalence 
between the matrix norms induced by these vector norms. 

Structured vector norm. We define a struc tured vector 
norm, following the approach presented in iRusso et’aH 
(l2013h . Specifically, assume x = [xi, • • • ,Xn\'‘ G 77” can 
be partitioned into K vectors x^ G 77”'“, k = I,-- - ,K, 

such that X = [x^ , • • • , with 

\x^\sk be the norm in K."''. Then, the structured norm of x 

is denoted by | • jc and defined as: 


NoteA„, = 


, which is un- 


-6.2859 2.066r 

^ 0.3562 0.1182y 
stable as it has eigenvalues of —6.3988 and 0.2311. Hence, 


_ Mi)+MA _ 
2 


\x\g = 

where the norm | • |s is defined in . 


( 21 ) 
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Given the same partition of x^, k = I,-- - ,K, consider 
another structured norm | • |g' based on using the vector 
norms | • |s/_ in k = 1, • • ■ ,K, and | • jg/ in such 


that \x\g' 




S'' 


induced by the vector norm | • |g in R” and ^s,r be that 
of I • Is ,, in R^. Then, it can be derived that fiG,ri^(j')) < 
^S,riA{r)) (IRusso et all l2013l) . Proposition 3 implies that 
P '12 P '21 can be the transactions coefficients between 

I • |g,i and | • |g .2 as well. 


Proposition 3 Let | • |g and | ■ |g' ke two structured norms 
defined as mentioned above. Let ts be the transaction coef¬ 
ficient from the norm | • |s' fo I ‘ |s ke the the trans¬ 

action coefficient from the norm | • |s' to | • |sj,. Then, we 

have \x\g' < Ts|17|s|a;|G with U = diag[Tk\^^i. 

This resu lt can be der i ved a s a consequence of those pre¬ 
sented bv iRusso et al.l (l2013h . 

IRusso et'aP (l2013h showed that contraction analysis can be 
used to carry out the hierarchical analysis and design of net¬ 
worked systems. Analogously, consider system (18). Parti¬ 
tion X into several sub-vectors, K vectors; x^ £ k = 
1, • • • ,K, so that X = • • • , x^^]^ with J2k=i = 

n. Each x^ corresponds to the k-th subsystem. Then, system 
(18) can be equivalently written in the following form: 


Therefore, suppose that switched system (18) is in the 
block-wise form (22). Then, Corollary 1, specifically in¬ 
equality (19), can be fulfilled assuming that it holds that 


(1/2) ^s,i(A(l)) -b /rs, 2 (^( 2 )) -b fr ■ (log/3(2 + log/32i) 


< 


—c. 


5 Synchronization in switched networks 

Finally, we consi der a network example in spired from one 
first presented in Idi Bernardo et ^ (l2014h . We assume the 
network equation is given by 


if = f{x''{t)) — kait)'^^^ Li jTx^ {t), z = (23) 

1=1 


x'== + = ,iT, (22) 

fc'=i 

where 


'^Aii(r) Ai2(r) • 

■ AiK{r)'^ 

^21 (x) A22 (x) • 

■ ^2Jf(x) 

\^ivi(x) AK2{r) ■ 

■ Akk^t)/ 


Here, x* £ M" stands for the state vector at node i, /(•) : 
R" —>■ R" the node dynamics, k is the coupling strength, 
L = [LifiY'j^i is the Laplacian matrix associated with a 
graph G = [V, E], where V = {1, • • • , m} is the node set 
and E the link set, by the way that for each Lij takes 
value — 1 if there is a link from node j to i and 0 otherwise, 
and La = — ^ij- ®’(^) takes values 0 or 1, implying 

that the diffusive coupling among the nodes in the graph G 
is only active when a{t) = 1 while it is not present when 
a{t) = 0. r £ R”’" stands for the inner coupling matrix. 
We can rewrite (23) in compact form as; 


and [Bj,--- ,B^]^ = B. 


X = E{x{t)) — k [<T(t)L 0 r] x{t) (24) 


The norm of x is defined by (21). Let Aijif) = \Aij{t)\ij, 
where the norm | • | y is defined by 

\Aij{r)\ij= sup \A,j{r)x^\s,. 

and consider the reduced K x K matrix 


i(r) = 


'^Aii(r) Ai2(r') 
A2i{r) A22{r) 


AiK{t) \ 
A2Kir) 


r = 1,2. 


V^ivi(x) AK 2 {r) ■■■ AKK{r)J 


Let I • |5_i and | • |5^2 be two norms in R*^ and P '12 and P '21 
be their transaction coefficients such that | • |5,1</312|-|S,2 
and 1-15.2 < P21 1 • Is.i- We take the multiple norms in R" as 
I • Ig.t = |^(x)|5,r, x = 1, 2. Let gG,r be the matrix measure 


wherex = Jx^^,--- £ R™”, 

E{x) = [/^(x^),--- ,/^(x"*)]^ and 0 is the Kronecker 
product. Assume that L is diagonalisable, i.e., there exists 
a nonsingular Q £ R™-™ such that L = Q~^JQ with a 
diagonal matrix J = where Aj, j = 1, • • • , m, 

are the eigenvalues of J, which are assumed to be real. 
Without loss of generality, we can assume that 0 = Ai > 
A2 > • • • > Am. Note that Ai = 0 is associated with the 
synchronization eigenvector [1, • • • ,1]^. 

Follow in g similar arguments to those by ICMrolL&J^gcora 
( 1991 ): Russo & di Bernard (1201 ih : Idi Bernardo et al. 
( 2014 ): lYi et al.l (12013h . synchronization of network (24) 
can be achieved if the following linear systems (obtained 
via linearisation and block diagonalization of (24)); 

p = [Df{w(t)) - kXi(7{t)]4>, i = 2, • • • ,m, (25) 

are contracting, where w{t) is a solution of the uncoupled 
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system 


w = f{w), (26) 

and D f{-) is the Jacobian of /(•). Assume 

T-L^: System (26) ha s an asymptotically stable attractor A (see 
lYi et al.Ll2013L Assumption 2 ). 



Definition 4 System (23) is said to synchronize if there 
exists (5 > 0 such that for any € iB(A,S)’ 

liiTLt^oo |x®(f) — x^(t)\ = 0 for all = I,--- ,m. 
Here, B(^A,S) = {y : inf^g^ jy — z\ < <5} denotes S- 
neighbourhood of set A 

According to Corollary 1, in particular, equation (16), we 
have the following result: 

Proposition 4 Consider two vector norms \ ■ |o and | ■ |i, 
which induce two matrix measures fJ-oi’) tmd /ii(-) respec¬ 
tively. Suppose thatHi holds and /ii(—F) < 0. If there exist 
To > 0 and c > 0 such that 


Fig. 3. The underlying graph topology. 



Fig. 4. Synchronisation dynamics of the network system (24) 
with random initial values that are picked from [—1,1]^, 

i = 1, 2, ■ • ■ , 10. Here, Err(t) = ^ 
with Xj = (1/m) 

To illustrate this result, we take 


- ^Mg{Df{w))%{nT,{nAl)T) 

+pii{Df{w) - kX 2 T)Ti{nT, (n + 1)T) 
+Mgi{nT, In + 1)T) log/3oi 


+^01 (nT, (n + 1)T) log/3io 


< —c 


(27) 


holds for all T > Tq and w G A where Tuis,t) stands 
for the duration of a(t) = u in the time interval (s, t] and 
Xfuv{s, t) is the number of switches from a{t) = uto a{t) = 
V, for all u V, u,v = 0A< then system (24) synchronises. 


Proof Note that for each i = 2, - ■ ■ , m, 
di{Df{w) - kXiV) 

< p.i{Df(w) - kX2T) + k{Xi - A2)/ri(-r) 

< di{Df(w) - kX 2 T). 


Due to ^i(—r) < 0, then (27) holds for alH = 2, • • • , m. It 
can then be proved in a straightforward manner that the vari¬ 
ational equations (25) are asymptotically stable for alH > 2 
by verifying condition (16) and following the same step s 
of the proof of Theorem 5.1 in Idi Bernardo et al.l (l2014h . 
Thus, followin g the a rguments in the proof of Theorem 
17 in lYi et al.l (1201 3l) . there is an invariant open convex 
set U G U D S, where 5 = {x = 

• • • , : x"^ = x^ G A,\/ i, j}, is invariant for the 

coupled system (24). Furthermore, S is an asymptotically 
stable set for system (24). This proves synchronization. ■ 


{ p{G[-wi + W 2 ] - g{wi)} 
G[wi - W 2 ] + UI 3 
-qw2, 


where g{wi) = mgWi-\-I/2(rm — mg){\wi-\-I\ — |mi —1|) 
as the Chua’ circuit with mg = —0.5, mi = —0.8, G = 
0.7, p = 9 and q = 7, which was reported to exhibit a 
double-scroll chaotic at tractor for the node dynamics (26) 
iMatsumoto et al.l(ll985h . Moreover, we take F as the identity 
matrix, and assume the underlining graph of 10 nodes has the 
structure shown in Fig. 3, associated with a Laplacian L with 
X 2 {L) = 2.7142. To study convergence, we choose weighted 
1-norms: |y|o = |?/i| + r 2 \y 2 \ + rglys], with r 2 = 3.4042 
and rg = 1.0369, and |y|i = y/|yip + |j/ 2 p + lyap, which 
implies (i) /3oi = 4.3163 and Ao = 1; (ii) do{Df (w)) < 
3.2829 for all t. We take fc = 1, so that y,i{Df(w) — 
kX 2 {L)r) < —7.4714, for all t. The switching signal a{t) 
is taken as follows: 

, , Jo tG [kT, kT + 1/4T) 

’ \l tG {kT+l/4T,{k + l)T) 

for some T > 0. Condition (27) in Proposition 4 implies that 
if T > 13.08, then system (24) synchronises. To illustrate 
this result, we take T = 14. Fig. 4 shows that, as expected, 
all nodes synchronise. 

6 Conclusions 

We have presented an extension of contraction analysis to 
switched Caratheodory systems. The key step was the def- 
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inition of switched matrix measures induced by multiple 
norms. Using these measures, it was possible to derive dif¬ 
ferent sets of sufficient conditions for asymptotic incremen¬ 
tal stability of the systems of interest. Most notably, it was 
possible to prove contraction and hence incremental stabil¬ 
ity by using different norms, each associated to a different 
mode of the switched system under investigation. This com¬ 
plements and extends in a highly nontrivial manner to the 
case of multiple nor ms, previous results presented by some 
of th e authors (see iRusso et"^ 1201 3t Idi Bernardo et al.[ 
l2014ll . where contraction was studied by using a common 
metric. The theoretical results were illustrated on a set of 
representative examples and applications showing the effec¬ 
tiveness of the proposed method. 
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